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ABSTRACT 

Noil-intrusive measurement* of scalar and momentum trans- 
port in turbulent wall plumes, using a combined technique of laser 
Doppler anemometry and laser-induced fluorescence, has shown 
some interesting features not present in the free jets or plumes. 
First, buoyancy-generation of turbulence is shown to be important 
throughout the flow field. Combined with low- Reynolds- number 
turbulence and near-wall effect, this may raise the anisotropic tur- 
bulence structure beyond the prediction of eddy-viscosity models. 
Second, the transverse scalar flaxes do not correspond only to the 
mean scalar gradients, as would be expected from gradient-diffusion 
modeling. Third, higher-order velocity- scalar correlations which 
describe turbulent transport phenomena could not be predicted us- 
ing simple turbulence models. 

A second-order closure simulation of turbulent adiabatic wall 
plumes, taking into account of the recent progress in scalar trans- 
port, near-wall effect and buoyancy, is reported in the current study 
to compare with the non-intrusive measurements. In spite of the 
small velocity scale of the wall plumes, the results showed that 
low- Reynolds- number correction is not critically important to pre- 
dict the adiabatic cases tested and cannot be applied beyond the 
maximum velocity location. The mean and turbulent velocity pro- 
file* are very closely predicted by the second-order closure models. 
But the scalar field is less satisfactory, with the scalar fluctuation 
level underpredicted. Strong intermittency of the low-Reynolds- 
number flow field is suspected of these discrepancies. The trends in 
second- and third-order velocity-scalar correlations, which describe 
turbulent transport phenomena, are also predicted in general, with 
the cross-streamwise correlations better than the streamwise one. 
Buoyancy terms modeling the pressure-correlation are shown to im- 
prove the prediction slightly. The effects of equilibrium time-scale 
ritio and boundary condition are also discussed. 


• Assist. Professor, ASME member 


INTRODUCTION 

Wall plumes are encountered above fires along surfaces, near 
baseboard heaters and electronic circuit boards, and in other con- 
fined natural convection processes. Turbulent wall plumes, al- 
though simple in geometry, are also a fundamentally interesting yet 
complicated flow. Wall jet has the characteristics of both boundary 
layers and free jets. Wall plumes are further complicated by buoy- 
ancy. It is characterized by intermittent, low- Reynolds- number and 
large-scale turbulence, involving the interactions of buoyancy, tur- 
bulence, and the damping effect of a solid surface. Therefore, pre- 
dicting this flow poses a great challenge. In addition, wall plumes 
provide an excellent opportunity for studying these interactions, 
since they are relatively thick, reducing problems of spatial res- 
olution, but do not exhibit the effects of large-scale disturbances 
(flapping) encountered in free plumes, due to the stabilizing effect 
of the wall. 

Past measurements of wall plume turbulence properties are very 
limited. Grella and Faeth (1975) used hot-wire anemometry to mea- 
sure streamwise velocity fluctuations in turbulent adiabatic wall 
plumes. Liburdy and Faeth (1978) and Liburdy et ah (1979) stud- 
ied turbulent thermal plumes along a vertical isotherm ad wall. Hot- 
wire anemometry was also used to measure mean and fluctuating 
velocities and temperatures and their correlations. However, these 
studies have several limitations, and no data are available to con- 
firm them. First of all, parasitic heat losses from thermal wall 
plumes are difficult to control; therefore, adiabatic conditions were 
only approximated by Grella and Faeth (1975). Furthermore, in- 
trusive probes are not very reliable for measuring the properties of 
wall plumes. Probes have large uncertainties when turbulence in- 
tensities are high, and near surfaces; therefore, properties are only 
measured reliably in a narrow region of the flow. Probes also dis- 
turb the flow, particularly near the edge, where flow reversals occur 
(Lai et al., 1986). 

Lai et al. (1986, 1987) reported the first non-intrusive mea- 
surements of mean and fluctuating velocities and scalars, as well 
as their correlations in turbulent adiabatic wall plumes. Problems 
associated with past experiments were eliminated entirely. Since 



the data of this experimental program is used to validate the SOC 
model in the present study, a brief summary is provided in the 
following section. These flows were analyzed using both mixing- 
length model and simplified k-t-f model with wall function but 
ignoring buoyancy/turbulence interactions entirely, similar to past 
analyses of buoyant round diffusion flames (Jeng et al., 1982; Jeng 
and Faeth, 1984). Both methods yielded satisfactory predictions 
of mean quantities in spite of the effects of low-Rc and large-scale 
coherent structures seen in the flow visualizations. However, tur- 
bulent fluctuations were underestimated. This was attributed to 
the effects of buoyancy /turbulence interactions by analysis of the 
production budget in the fc-eqaation using velocity /scalar correla- 
tions. Due to the small velocity scale of wail plumes, the y + at the 
maximum velocity point y m is well inside the buffer layer (ca. 18); 
even at the outer edge of the flow y + is very small (ca. 150, using 
u y as defined by Kraka and Eskinari (1964) for wall jet). This puts 
the first grid literally within the laminar sublayer. Therefore, in 
spite of the apparent success, implementation of the wall function 
in k-of model for wall plumes is debatable. 

However, interesting features of the flow properties remain be- 
yond the prediction of the eddy- viscosity models, such as mixing- 
length or k-€-f models. For example, the anisotropic turbulence 
intensities in the stream wise (open symbols) and transverse (closed 
symbols) direction cannot be predicted. Even no counter- gradient 
diffusion is observed in the measurements, the transverse scalar 
fluxes do not correspond to the mean scalar gradients, as would be 
expected from the gradient- diffusion assumption in eddy- viscosity 
models (Lai, 1985). Therefore, the diffusion models based on the 
first principle as suggested by Lumley for SOC model is preferred 
over the simple gradient- diffusion assumption for the eddy- viscosity 
models. 

In addition, higher-order velocity-scalar correlations which de- 
scribe turbulent transport phenomena could not be predicted using 
simple turbulence models. Triple moments dominate the diffusion 
terms in the second moment equations; for example, 

</tMtj> and <pu±> appear in the second- moment equations for 
the Reynolds stress <UtU ; >, scalar flux </t^> and scalar variance 
</*>. These correlation measurements would therefore validate 
and provide new insight into modeling of the scalar transport pro- 
cesses. 

One possibility is to apply the algebraic stress model, which is 
considered a modified version of the i-c model without the gradient- 
diffusion assumption (Rodi, 1972, and Gibson and Launder, 1978). 
Therefore, the algebraic model could have a wider range of appli- 
cability than the eddy- viscosity models, although the assumption 
used in deriving this model needs further justification. There are al- 
ready some buoyant flows considered with the algebraic model. For 
example: buoyancy- modified free surface flows (Gibson and Laun- 
der, 1976, Hossain and Rodi, 1977, and McGuirk and Papadim- 
itriou, 1985); natural convection on a flat plate and a heated cavity 
(Humphrey and To, 1985, 1986, and To and Humphrey, 1986); and 
a developing buoyant plume (McGuirk and Papadimitriou, 1985). 
However, no complete second-order closure (SOC) simulation has 
been performed, on the two-dimensional wall plumes. 

The recent renewed effort in SOC models, aiming at improv- 
ing the physical ingredients of second-moment closure schemes, as 
reviewed by Launder (1989), holds promise for improving the re- 
liability with which we can predict industrially interesting flows. 
Unlike previous models, second-order closure models directly solve 
the transport equations for Reynolds stresses and scalar fluxes. In 
order to close the system equation, a number of terms must be mod- 
eled, such as the third order moments, the pressure- correlations, 
and the transport terms for dissipation rates. Various models have 
been proposed for these terms (Daly and Harlow 1970, Hanjaiic 


and Launder, 1972, and Lumley, 1978). The first principle- based 
diffusion terms by Lumley has been shown to produce more realistic 
behaviors in wall plumes (Lai and Jeng, 1990). Current develop- 
ment is primarily attributed to Lumiey’s extensive article (1978), 
which is based on tensorial algebra, realizability concept and first 
principle modeling of diffusion terms. Further impetus is given by 
the direct numerical simulation efforts that have been emerging 
over the past five years from the Stanford /NAS A Ames Research 
conglomerate. Noteworthy developments are in progress contribut- 
ing to the area of scalar transport, such as: pressure-correlation, 
diffusive transport, dissipation and near- wall effect. 

The whole area of low- Re near- wall turbulence is one where 
the scalar flux development lags well behind that of the Reynolds 
stresses, even for flow parallel to a wall. A consistent low- Re treat- 
ment of both the velocities and scalar is currently in progress (Shih 
and Mansour, 1990, 1992, Shima, 1988, and Launder and Shima, 
1990), aided by direct numerical simulation data (Kim and Moin, 
1987, and Rogers et al., 1986). The goal of the present study is, 
therefore, to validate the performance of SOC models on a turbulent 
adiabatic wall plume, based on Lai’s (1985) data. The model con- 
sidered incorporates recent developments in modeling scalar trans- 
port, pressure- correlation, low- Re and buoyant flows. 


THE EXPERIMENTAL PROGRAM 

This section briefly summarized the experimental program, with 
emphasis on the combined velocity- scalar measurements. A full ac- 
count of the complete experimental program is provided by Lai 
(1985). The test flow was generated by weQ-mixed carbon diox- 
ide/air flowing downward along a plane wall in still air. The carbon 
dioxide/air mixture entered a plenum at the top of the wall and 
then passed through a series of baffles to achieve a uniform two- 
dimensional flow at the 21-mm-high exit slot. The test wall was 
1000mm long and 800mm wide, and had 305- mm- high side walls 
to help preserve two-dimensionality. Windows in the side walls pro- 
vided optical access for structure measurements at zjb = 0 1 (the 
initial condition), 10, 20, and 37.5. 

Gas flow to the top of the wall was provided by an oil free air 
compressor and commercial-grade carbon dioxide stored in cylin- 
ders. Flows were controlled by pressure regulators and metered 
by critical-flow orifices. After mixing, the flow passed through a 
25mm(id.) x 37m long tube to ensure complete mixing and equi- 
libration to local room temperatures. Uncertainties in flow mea- 
surements and initial gas composition were less than three percent. 
In the experiments reported here the density ratios of the carbon 
dioxide/air mixture (at the exit slot) to the ambient air were 1.02 
and 1.04. The mean velocities at the exit slot were set at 0.31 and 
0.43m/s respectively, to match an asymptotic Froude number of 5, 
which is similar to the one measured by Greba and Faeth (1975) for 
vertical adiabatic wall plumes. This is confirmed by the fact that 
the maximum mean velocities measured at the furthest downstream 
position are almost the same as those at the slot exit. Therefore, 
the flow neither under-nor over- accelerated at the slot exit, which 
reduces the distance for flow development. 

A combined laser Doppler anemometry and laser-induced fluo- 
rescence (LDA/LIF) system was used to measure the velocity and 
scalar (mixture fraction) simultaneously; see Lai and Faeth (1987) 
for a detailed description of the arrangement. Iodine vapor was 
seeded in the flow to provide the LEF signal, while aluminum oxide 
particles (around 500 am diameter) were seeded in both the slot 
ajid the ambient gases for the LDA system. The LIF signal was 
separated from light scattered from the laser line using long- pass 
optical filters. The LDA and LIF were both sampled for 100 sec- 
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on ds at 80ff i using a low-pass anti-aliasing filter with a 40 Hz cutoff 
frequency. Measurements of power spectral densities of scalar fluc- 
tuations showed that power spectra were on the order of 0.1 percent 
of maximum value at 40#* at the worst (Lai, 1985); therefore, the 
frequency response of this system was adequate. The signals were 
processed to yield time averages; however, the distinction between 
FavTe and time-averaged quantities is small for these flows (less 
than five percent). 

The measurements are estimated to have the following experi- 
mental uncertainties (95 percent confidence): mean and fluctuating 
velocities less than four and six percent, mean and fluctuating con- 
centrations less than five and ten percent, Reynolds stress less than 
twenty percent, and turbulent mass fluxes (second-order correlation 
of velocity and scalar) less than fifteen percent. These estimates are 
based on the maximum value of each quantity and are proportion- 
ately higher elsewhere. Higher-order correlations have even larger 
uncertainties, which can be estimated from these estimates using 
conventional error- propagation analysis. 

The two-dimensionality of the plume was checked by computing 
momentum and buoyancy fluxes along the wall, si m i l a r to Launder 
and Rodi (1981). Mean momentum was conserved within ten per- 
cent and buoyancy flux was conserved within five percent (Lai, Jeng 
and Faeth, 1986), although due to turbulent fluxes in variable den- 
sity, the uncertainty of this conservation check is found to be about 
fix percent for momentum flux and four percent for buoyancy flux 
(Lai and Faeth, 1987). Aspect ratios, based on displacement, mo- 
mentum, scalar mixing, and analogous flow widths, were all greater 
than 10, also suggesting conditions conducive to two-dimensional 
flow. 

Initial profiles of stream wise and fluctuating velocities were mea- 
sured, and they were similar for all the flows. Ambient stratification 
could not be reduced below 0.26 °K/m without unduly disturbing 
the desired stagnant ambient environment of the wall. 


THEORETICAL METHODS 

The analysis treats a steady two-dimensional vertical wall plume 
in a stagnant environment having constant properties. Boundary 
layer approximations are applied, and dissipation and kinetic en- 
ergy are ignored in the governing equations for mean quantities, 
with little error. Carbon dioxide and buoyancy are conserved. Typ- 
ical of most analyses of turbulent mixing, the exchange coefficients 
are only a function of mixture fraction F (the fraction of mass which 
originated from the slot). Therefore, passive scalar properties are 
only a function of the mixture fraction, which is a conserved scalar 
of the flow. 

Turbulent properties were treated using two second-order clo- 
sure models. Here we consider the model representative of current 
development. The second-moment closure model is primarily based 
on the works by Lumley (1978); the scalar transport, Shih, Lumley 
and Chen (1990); the near- wall treatment, Shih and Mansour (1990, 
1992); the pressure- correlation terms; Shih, Shabbir and Lumley 
(1991). The low - Re SOC model is validated using direct numer- 
ical simulation of the fully -developed channel flow. No change of 
model coefficients is attempted to fit the predictions to the data. 
In the following section, only the results of model derivations are 
■ummarized; details of the derivations can be found in the cited 
references. 

Based on the above assumptions, a system of eleven transport 
equations (in addition to the conservation equations of mass) are 
to be solved. These include two conservation equations of mean 
properties (velocity and mixture fraction F, five modeled equations 
associated with turbulent momentum transport (Reynolds stress 


terms <u 3 >, <u 2 >,<w J >, <ut7> and <e>, which is the dissipa- 
tion rate of turbulence kinetic energy <!:> ), and four modeled 
equations associated with turbulent scalar transport (variance of 
mixture fraction < f 2 > , its dissipation rate <c^>, and scalar fluxes 
</u> and <fv>). 


Formulation 


Invoking the boundary- layer assumption, the mean continuity, 
momentum, and mixture fraction equations for the current wall 
plume calculations are: 


9pU dpV 
dx dy 
dpUU dpUV 
dx dy 
dpUF dpVF 
dx + dy 
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and the density is defined by the state relationship: 
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where g is the acceleration of gravity, po is the density of the fluid 
at the slot, and is the ambient density. Upper case and < > 
symbols stand for the ensemble- averaged properties and lower case 
symbols the fluctuating properties. 

For simplicity, tensorial notation is used for the following mod- 
eled equations. 


U k <UiUj> j = -<u,u ,->17** - <u 3 u k >Ui,k - [<«»*>**> 


- tij + <pu,>g x /p + <f*it>g 3 /p - ( B k} , + Bin^gkfp ( 5 ) 

U k <f*i>jt = -<fuj>Uij - <9i*j>Fj 

- {<fuiVj> - + 2U j, k T ijk 

- IIi<c>/<v s > + <pf>9i/P - B^gk/p (6) 

U k <f*>* = -2 <fu j >F J - </S > j ~ 2« ; > (7) 

U k <€>Jt = -<CU j>j -$<€><€>/ <q*> (8) 

u k = -<e'uj>j - ip 1 «*>*/</*> ( 9 ) 

where 

<^> = <u,u,> ( 10 ) 

<c> and <(J> are the dissipation rate of <q J > and <P> : 

«> = <e»>/2 = i/<u (11) 

<t*> = a<fjcfjc> (12) 

and 

(i3> 


c and c 7 are the model coefficients as 0.2 for the pressure stirring 
terms (Lumley, 1978). 

Attributing to Rotta's (1951) pioneering work, two kinds of in- 
teractions contribute to the pressure straining effect, one involving 
only fluctuating quantities and the other arising from the pres- 
ence of the mean rate of strain. These two iterations are termed 
the Return- to- Isotropy (or Slow) part and the Rapid part. In the 
above equations, the Rapid terms (/,)*/ and fi jk ), the Return-to- 
Isotropy terms (D^ and Hi), the third- moment terms (<u,Ujtijt>, 
</«,«)>, <P*i> and <*V>), the transport terms in the 

dissipation rate equations p and and the buoyancy terms doe 


3 


to pressure* correlation ( B tJ k and B{ } ) need to be modeled. The 
rapid term in the Reynolds stress equation U modeled by Sbih and 
Lumley (1990) and Shah et al. (1991), as 


<«*>, 




W = ( 4( » s r ~ e n ( <> ~ Mw) 3 


WpKj ~ S ?M 

- + t'jhpi + + Mu ~ yMw ~ 3 *w 4 «') 

+ - 34 w 6,7 - (14) 


where 


hij = <u,Uj>/<^> - $»i/3 

C.i = 1 + 0.8{1 - exp[(-fi t /40) 3 )}F/ /3 


c. 3 = i - f; 




a. = 


<«»>* 
9 v«> 


and 


F/ = 1 + 27 III + 9// 

m = 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 


+ C03 Z> + bjjfc </**«>) “b CDzbikb)cp<fUp> 


where 


The expression of R t in C Q \ is to account for the wall effect and the 
low- Re nature of the flow (Shih and Mansoor, 1990, 1992). 

The scalar counterpart of the rapid term is modeled as 

lijk = ^«> < / u k > " ^ <f u %>) + Cd\ bij<fvk> 

( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 


Cd\ = 

Cm = 

Cm- 7 + J 7 


and 


Fd - j — 2^ 3 ^ 3 + 

“ </ 2 ><^> - <fu r ><fup> 


Model constants cj, e 3 and Ci are taken to be 1.8, -1.8 and 4.5 
respectively (Sbih et al., 1991). 

The return-to- isotropy term is modeled together with dissipa- 
tion rate (Shih and Mansour, 1992), as 

A- 3 * - €ij = - «>0 “ /») " <(? 2 > [ 2<u,u ^ > 

+ 4(<u,ujt>n J ni + <Ujtifc>n,rn) + 2<u*Uf>n*n/n t n ; ] (28) 

where 


fi = 2+f exp(-^)[^ + 80.1 ln[ 1 + 62.4(-// + 2.3///)] j 


K“) l*, 1 ' 

/« = exp | 


(29) 

(30) 


and 

Rr = UrS/v (31) 

where U r is the shear velocity, 6 is the characteristic length scale 
taken to be distance from the wall to the maximum velocity point 

*»■ 

The scalar counterpart of the slow term is modeled as 


where 


*/ = £ + r- 


n. = <t/<jui> 


(ft - 2)(1 — d,,<f,j)/12 


yj — Jji _ . a I 

(1 — )/6 + d<j bjkdkt — b > jd%j 


n 


(32) 


(33) 


ud 


// = 1.1 + 0.55{/J - 1) tanh[4(r - 1)] (34) 

r _<jW> , 

<e> <}> 

The right- hand -side of the dissipation rate equations is modeled as: 


<€Uk> = - 


5(40 


(36) 
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2.4 

-0.15(1 - Fj) 

3.8 
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(38) 

(39) 

[1 - 0.331n(l -55//)] (40) 

(41) 

(42) 

(43) 

(44) 


*-?GH 

V)(t)“ 

(2+ -) |l - 0.1 (l - y)j for r>r, (45) 

where the eqofljbriam time-scale ratio r t is taken to be 1.6 
(Shih, 1990). These modeled terms are found to yield much better 
results compared with those using the model proposed by Hanjaiic 
and Launder (1976). 

The third-moment terms modeled by Lumley (1978) are used 
in this study. 

<UiUjU k > = ~ 3fl <t> ( <tttU >’ ><u,u,> J 1 


+ + <«,Up><u J u t > J> ) 




9/J 


(46) 


4 


(47) 

</«.«;> = ~ Jjj + ^) <(> { ' <U> ' ii><! ' l ' >k + <w></u '^ 

+<n,Uj> k <fvk>) + 3 (^+ 2 ^) <?3/> (48) 

<ff> = " (2+^fe> (2<U, “ ,></U7> * + <,?7> » </u ’ >) 

(49) 

</,u * >= “( 57 wk^ {2<!Kr><fUk> * + < ' i * ><p> * ) 

(50) 

Tie simpler formulations of Hanjalic and Launder (1976) for 
<UiUjUk> and <tuj r > were also tested; the results are very similar. 

Finally the buoyancy terms Bij^ and Bf } modeled by Shih et al. 
(1991) are used in this study without alteration. Their expressions 
ran be found in the reference. 


Computation 

A parabolic computer program similar to GENMIX (Spalding, 
1980) was used to solve the highly- coupled parabolic partial differ- 
ential equations. However, this program solves modeled transport 
equations simultaneously rather than sequentially, by using a block- 
tridiagonal solver. This is to alleviate the problem of numerical 
instabilities due to the complex coupling of the modeled equations 
associated with second-order closures. The numerical schemes use 
a non-staggered grid system to facilitate wall boundary condition 
treatment. A semi- Implicit scheme is applied as the integration is 
performed downstream. The scheme has a second-order accuracy in 
the transverse direction and a first-order accuracy in the marching 
direction. 

The boundary conditions are all set identically to zero except 
for F and <e> at the wall: 

U = <u 3 > = <t> 3 > 3= <w*> - <uv> - 0 


<p> = </«> = <fv> = <€*> 


dF 

dy 


= 0 ; <€> = 


v ( d<q 3 

4<g 3 > \ dy 




= 0 

at y - 0 


(51) 


For the results presented, fifty grids point was used inside the plume 
layer. The measured conditions at the exit slot (Lai, 1985) were 
used as initial conditions in the computation. 


RESULTS AND DISCUSSIONS 

Since the data for wall plumes are similar, all data used for 
comparison is for the test case p 0 /p Q © = 1 -02. Only the data at 
the furthest downstream position z/b = 37.5 is presented. For con- 
venience, all plotted data and predictions are non- dimensional) zed 
with respect to the maximum mean properties (U m and F m ) at 
the particular downstream location, x/b = 37.5 . To demonstrate 
the effect of the low- Tie treatment, three different implementations 
were used for comparison. First one is the complete low-Re SOC 
model as formulated above except for the ^ in Eq. (41) is change 
from 2.4 to 2.1, as suggested in Shih and Man sour (1990); this is 
designated as the LL model. Second one, which is designated as 
the HH model, does not use low- Re treatment at all; i.e., R< in Eqs. 
(16) and (29) is set to be infinity and in Eq. (42) is set to be 
zero for the entire region. Third one, the LH model, implements 


low- Re treatment only ap to y m , beyond which the regular high Re 
model is used. 

Figure 1 shows the comparison between the predictions and 
measurements for the velocity structure. Interestingly, the HH and 
LH model predictions axe very similar, with the LH slightly better. 
Both the HH and LH results agree very well with the experimental 
data, only underpredicting the plume width slightly. This may be 
why that Lai and Jeng (1990), using simple boundary conditions 
In earlier SOC model (Lumley, 1978) can predict wall plume with 
relative success. Both models predict the correct maximum turbu- 
lence intensities and Reynolds stress, which is & major success not 
achieved in previous analysis (Lai et al., 1985-1987). The LL model 
results do not agree well with the experimental data It overpre- 
dicts the plume width on the free side, while underpredicting the 
turbulent intensities, suggesting that the turbulence is overdajnped . 
Therefore, Low- Re treatment should not be applied beyond y m for 
wall plumes. 

Figure 2 shows the comparison between the predicted and mea- 
sured scalar structure. In general, the agreement is not as well as 
the velocity structure. The change of sign in the curvatures of F 
around y m was more pronouncedly predicted than the original data 
as reported by Lai (1985). The trends in second- and third-order 
velocity-scalar correlations, which describe turbulent transport phe- 
nomena, are also predicted in general, taking into the account of 
the larger experimental uncertainty of these measurements. How- 
ever, the predicted scalar fluctuations are notably lower than the 
measured levels, with the stream wise correlations worse than the 
cross-streamwise ones. Underpredicting < / u> will underestimate 
buoyancy production (Lai and Faeth, 1987) and therefore under- 
estimate the plume width. When comparing the velocity- scalar 
correlation coefficient results, (i.e., when normalized with respect 
to the local turbulence quantities <tq> and </>, rather than U m 
and F m the agreement is much better for all the correlations in 
the central part of the flow. However, small numerical values of 
the predicted fluctuations at the outer edge will sometimes cause 
the correlation coefficients to overshoot beyond unity and violate 
the Schwarz Inequality. This suggests that the combination of the 
submodels in the SOC is not completely realizable at the edge; how- 
ever, this should not cause the laxge discrepancy observed in the 
scalar results. The nature of the LDA/LIF measurement may have 
introduced slightly larger fluctuations in the scalar measurement, 
since the scattered light from LDA seeding particles may stiff fil- 
ter through the long-pass filters and contribute to the fluorescence 
signal. However, this error has been shown to be small (Lai and 
Faeth, 1986). A more plausible explanation can be found in the flow 
visualizations picture and the probability density function (pdf) of 
the scalar measurements (Lai, 1985). They show large structure of 
relatively unmixed eddies and deep intermittency spikes penetrat- 
ing very close to the wall. This may be also why the Low- Reynolds 
number treatment is not critically important as shown previously. 

For comparison, low-Re fr-e- / predictions of the wall plumes 
are also carried out. If low-Re correction is used instead of the waff 
function for the Jfc-e-/ model, one can integrate directly to the waff. 
The first version of the low-Re k-( mode! was proposed by Jones and 
Launder (1972) to predict relaminarization in severely accelerated 
boundary layer. Since then several authors have proposed alternate 
versions of low-Re h-t models. Patel et al. (1985) and Bernard 
(1986) reviewed the performance of the these models for boundaxy 
layer flows and found comparative results. All these models have 
the following basic form: 

U,kj = [(*' + k j\ + + U„)-<-D (52) 


5 


CONCLUSIONS 


Uj'j = [(* + tj] + Ca/i^UijiVij + 

-C <2 hj + E (53) 

where 

* = (54) 

The model of Chien (1982) is selected for the present study, since it 
contains only simple algebraic functions and do not involve differ- 
ential operators and therefore easier to implement. The constants 
and functions for this model axe: 

= 0.09; C < i = 1.35; C a = 1-8; a k = 1.0; a t - 1.3 (55) 

= 1 - exp(-0.01 15y + ) (56) 

fi = 10; h = 1 - 0.22exp[-(7?i/6) 3 ] (57) 

D = 2i/fc/y 3 ; E = -2K«/y 3 )«p(-0.5y + ) (58) 

For simplicity, buoyancy /turbulence interactions were not consid- 
ered; only mean momentum equation include the body force term. 
The computation is performed with fifty grid points in the trans- 
verse direction. Turbulent scalar transport term is closed by relat- 
ing to the turbulent viscosity through a constant turbulent Prandtl/ 
Schmidt number of 0.9. For comparison, implementation of the wall 
function is also performed by setting u + = y+ since y + is very small 
near the wall. Figure 3 shows the comparison of measurements and 
predictions by k-c-/ models. The result shows that the low- .ft* 
model overdamps the plume development beyond y m , while the 
results using wall function is better outside, but the development 
inside y m is underdamped. The predicted mean mixture fraction 
it a simple profile as suggested by the gradient-diffusion assump- 
tion. The implementation of buoyancy term can not be exact as in 
the SOC. A switch of the Low- .Re treatment in k-«-/ near will 

cause the computation to be unstable. These demonstrate that the 
limitation of eddy- viscosity isotropic models to predict the flows. 
However, the performance of turbulence models which are generally 
validated with high- Reynolds number flows may be limited in their 
performance in wall plumes. 

Figure 4 shows the triple correlation of velocities for the three 
SOC models. The trend and level of these correlations are similar 
to those found by Dekeyster and Launder (1983). The profile for 
the LL model is notably smaller than the other two due to damp- 
ing of the wall. Without the Low-fte treatment, <uuv> is larger 
<wv> for the HH model and part of the LH model results. Fig- 
ure 5 plots the time-scale ratio profiles. Both LL and HH model 
show a large variations in the time scale ratio except the LH model, 
suggesting the deficiency of the LL and HH models in predicting 
this flow. The equilibrium time-scale ratio r e is difficult to define 
theoretically, and the form Eqs. (44) and (45) of is derived crudely. 
Therefore r e is considered to be an adjustable constant. Changing 
the r e from 1.6 to 4.0 for the LH model shows a more even curve; 
however, the improvement in the prediction is only slight. The ef- 
fect of the pressure- correlation buoyancy term in shown in Figure 
6. Excluding this extraneou&iy modeled term deteriorate the per- 
formance slightly. However, dropping the buoyancy term due to 
velocity-scalar correlation, which is exact, will degrade the predic- 
tion significantly. The measured non- zero boundary conditions for 
<u>, <t;> and </> suggest that a zero-gradient boundary condi- 
tion for these properties could improve the prediction, as shown in 
Figure 7. 


A second-order closure simulation of turbulent adiabatic wall 
plumes, taking into account of the recent progress in scalar trans- 
port, near-wall effect and buoyancy, is reported in the current study 
to compare with the non-intmsive measurements. In spite of the 
small velocity scale of the wall plumes, the results showed that 
low- Reynolds- number correction is not critically important for the 
second-order closure model to predict the adiabatic cases tested 
and cannot be applied beyond the maximum velocity location. The 
mean and turbulent velocity profiles are very closely predicted by 
the second-order doeure models. But the scalar field is less sat- 
isfactory, with the scalar fluctuation level underpredicted. Strong 
intermittency of the low- Reynolds- number flow field is suspected of 
these discrepandes. The trends in secoad- and third-order vdocity- 
scalar correlations, which describe turbulent transport phenomena, 
are also predicted in general, with the cross- stream wise correlations 
better than the streamwise one. Buoyancy terms modeling the 
pressure- correlation axe shown to improve the prediction slightly. 
The equilibrium time-scale ratio for the wall plumes is found to be 
around 4. Zero-gradient boundary conditions for the fluctuating 
velocities and scalar also improve the prediction slightly. 
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NOMENCLATURE 


Biik 

BU 


b 

K 

c, </ 

Cji i 

Cl,C 2 ,C 3 

D 

E 

F 

Fd 

Ft 

f 

U 

fuh 

L 

9* 

hjki 

li,k 

II 

III 
k 

7 2 

Re 

Rt 

Rr 


buoyancy term due to pressure correlation, Eq.(5) 
scalar buoyancy term due to pressure correlation, 
Eq.(6) 

exit dot height 

Eq.(15) 

model constants, in Eq.(5) and (6) 

model coeffidents, in Eq.(23)-(25) 

model constants, in Eq.(61) 

modd constants, in Eq.(23)-(25) 

modd coefficients, in Eq.(16) and (17) 

near-wall treatment in k-equation (58), £q.(63) 

near-wall treatment in Eq.(27) 

near-wall treatment in e-equation (59), Eq.(63) 

mean mixture fraction 

function, Eq.(26) 

function of invarients, Eq.(19) 

fluctuating mixture fraction 

damping function for the slow term, Eq,(30) 

modd coefficients, Eq.(62) 

damping funtion for Low- ft* k-f modd, Eq.(62) 

acceleration of gravity 

Rapid term, Eq.(14) 

Scalar Rapid term, Eq.(22) 

In varient, Eq.(20) 

In varient, Eq.(21) 
turbulent kinetic energy 

/ovuj or twice the turbulent kinetic energy, 2 k 
Reynolds number 
Reynolds number, Eq.(18) 

Reynolds number, Eq.(31) 
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r time scale ratio, Eq.(35) 

r € equilibrium time scale ratio, Eqs.(44) and (45) 

[J , V mean velocity in (z,y) direction 

U % mean velocity in the i-direction 

U T shear velocity 

fluctuating velocity in the i-direction 
u, v, to fluctuating velocity in the (z, y, z) direction 

x downstream distance from the slot exit 

y normal distance from the wall 

y + nondimentional normal distance from the wall, 

Vry/» 

y_ normal distance from the wall with maximum U 


0 

e 

<•> 


Vt 



p 

Po 

0k> 

¥ 

1 ^ 0 , 01 , 02,^3 

0o» 0i 
¥ 


diffusion coefficient 
Eq.(29) 

Kronecker delta 

dissipation rate of turbulent kinetic energy, k 

dissipation rate of mixture fraction variance, / 3 

dissipation rate tensor of Reynolds stress, 

dynamic viscosity 

kinematic viscosity 

turbulent kinematic viscosity, Eq.(60) 

slow term, Eq.(28) 

scalar slow term, Eq.(32) 

density of the fluid 

density of the fluid at the slot exit 

turbulent P ran tl /Schmidt numbers, Eq.(61) 

Eq.(35) 

Eq.(38) 

model coefficients Eq.(40)-(43) 
model coefficients Eq.(44)-{45) 

Eq.(39) 


Subscript 



maximum mean properties 
ambient condition 


«n 

dXi 


Superscript 

( * ) ensemble average 

<> ensemble average 
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